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Dedicated to Mike Eastwood on the occasion of his 60th birthday. 

Abstract. We discuss the twistor correspondence between path geometries in three di- 
mensions with vanishing Wilczynski invariants and anti-self-dual conformal structures of 
("■^ , signature (2, 2). We show how to reconstruct a system of ODEs with vanishing invariants 

rvQ ' for a given conformal structure, highlighting the Ricci-flat case in particular. Using this 

. ,^ , framework, we give a new derivation of the Wilczynski invariants for a system of ODEs 

flj . whose solution space is endowed with a conformal structure. We explain how to recon- 

struct the conformal structure directly from the integral curves, and present new examples 
of systems of ODEs with point symmetry algebra of dimension four and greater which give 
CO ' rise to anti-self-dual structures with conformal symmetry algebra of the same dimension. 

Some of these examples are (2, 2) analogues of plane wave space-times in General Rela- 
tivity. Finally we discuss a variational principle for twistor curves arising from the Finsler 
structures with scalar flag curvature. 
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r-| ' 1. Introduction 

-(— > \ 

^ A path geometry on an open set 7" C W" consists of unparametrised curves, one curve 

through each point in T in each direction. Regarding these curves as solutions to a system 
of (n — 1) second order ODEs, one can give an alternative definition of the path geometry as 
cn ! an equivalence class of systems of second order ODEs, where two systems are regarded as 

^ ' equivalent if they can be mapped into each other by a change of dependent and independent 

variables. In this paper we shall study three-dimensional path geometries encoded into a 
system of two second order ODEs 

^. Y" = F{X,Y,Z,Y',Z'), Z" = G{X,Y,Z,Y',Z'), (1.1) 

(^ . where Y' = dY/dX etc, and {F, G) are arbitrary functions on an open set in M^ which we 

^N I assume to be of class C^ . Two such systems are locally equivalent if they are related by a 

point transformation {X, Y, Z) — )• {X{X, Y, Z), Y{X, Y, Z), Z{X, Y,Z)). A natural question 
arises: Given a system (jl.ip is it equivalent to a pair of trivial ODEs Y" = 0, Z" = 0? 
K> I The answer comes down to constructing a set of invariants for (jl.ip . These invariants can 

^ ■ be divided into two groups C and V. Vanishing of the invariants from each of these groups 

selects a subclass of three-dimensional path geometries. 

• The system (jl.ip belongs to the conformal branch if its four-dimensional solution 
space M admits a conformal structure such that the two-dimensional surfaces in 
M corresponding to points in T are totally isotropic (this condition uniquely de- 
termines the conformal structure). The conformal branch is characterised by the 
vanishing of three fundamental Wilczynski invariants C given by the expressions 
()2.ip . Systems with vanishing Wilczynski invariants are also referred to as torsion- 
free [17j and we will use this terminology several times throughout. 

• The system (jl.ip belongs to the projective branch if its integral curves in T are 
unparametrised geodesies of some projective connection. The projective branch is 
characterised by the vanishing of eight invariants V given by the expressions ()2.2p . 
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A system is point equivalent to a trivial one if and only if it belongs to both the conformal 
and projective branches |17| . In this case the symmetry algebra of (|1.1|) is isomorphic to 
5[(4,R). The corresponding Lie group P5L(4, R) acts projectively on T = MP preserving 
the unparametrised geodesies of the flat projective connection. 

In this paper we shall concentrate on the conformal branch, and make use of the local 
isomorphism 

P5L(4,M)^50(3,3), 

where 5(9(3, 3) is the conformal group of the flat conformal structure on the solution space 
M. If the curvature of the conformal structure does not vanish, then the conformal sym- 
metry group is a proper subgroup of 50(3, 3), and conformal Killing vectors on M give rise 
to point symmetries of the system (jl.ip (see Lemma 16. ip . We shall exploit this correspon- 
dence to construct several examples of torsion-free ODE systems and the corresponding 
conformal structures which admit a large symmetry group, 9-dimensional symmetry being 
the sub-maximal case. 

In the rest of the paper we shall call T the twistor space, and the integral curves of (jl.ip 
the twistor lines. In the next section we shall give explicit expressions for point invariants 
of (jl.ip . In section [3] we summarise the twistor correspondence anti-self-dual conformal 
structures which underlies all constructions in this paper. In section HI we describe how, 
given an ASD conformal structure of (2, 2) signature, to construct a system of ODEs on 
the twistor space T whose integral curves are the twistor lines. Using this framework, in 
section [5] we present a new derivation for the condition of vanishing Wilczynski invariants 
assuming solely that the space of solutions be endowed with a conformal structure. We 
demonstrate how this construction proceeds in the Ricci-flat case, and show (Theorem 15. ip 
that in this case the system of ODEs can be read-off directly from the solution of Pleban- 
ski's heavenly equation. In section [6l we explore the local isomorphism between groups of 
point symmetries of (jl.ip and conformal symmetries of the corresponding conformal struc- 
ture. In section \7\ we construct examples of torsion-free ODE systems in three dimensions 
with point symmetry algebras of dimensions between nine and two. Finally in section [8] we 
explore the connection between systems of ODEs with vanishing Wilczynski invariants, and 
unparametrised geodesies of Finsler structures with scalar flag curvature. This correspon- 
dence gives rise to yet another class of examples. Some of these examples have gravitational 
analogues in the theory of plane wave space-times, and the whole construction sheds new 
light on variational aspects of the Nonlinear Graviton theorem. 

It is a pleasure to dedicate this paper to Mike Eastwood on the occasion of his 60th 
birthday. Mike has made important contributions to twistor theory and conformal geometry 
over the last thirty five years. He and his coworkers [3] have developed an approach to 
symmetry of differential equations which underlies some of our work. 

Acknowledgements. We are grateful to Boris Doubrov for telling us about the torsion-free 
class of systems of ODEs, and for many useful discussions. In particular example (j7.2p 
was provided by Boris. We also thank Thomas Mettler for pointing out the connection 
between the torsion-free systems and Finsler structures with scalar flag curvature, and 
thank anonymous referees for pointing out several inaccuracies in the original manuscript. 

2. Systems of ODEs and their invariants 

To present the invariants in a compact form set Y = {Y, Z) and F = {F,G), where 
the capital letter indices A,B,C, . . . take values 0, 1. Moreover we set 

,A _ dY^ _^ d _ d , . d , . d 



P =^^T^, and ^T7 = T7T7 + -P T^TTT + i^ 



dX ' dX dX dY^ dP^' 
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The three Wilczynski conditions C are given by 

T - lTrace(r)]I = 0, where T^ = -|^ - 1^|^ + 1 A|^. (2.1) 

2 ^ ^ ' ^ dY^ 4 dP^ dP^ 2dX dP^ ^ ' 

Our terminology is motivated by the work of Doubrov [8t i9j , who has demonstrated that the 
conditions ()2.ip for systems of non-hnear ODEs are equivalent to the classical Wilczynski 
invariants [30] on the linearisations of these ODEs. 

The complementary set of invariant conditions V characterising the projective branch is 

S (BCD) = (2.2) 

where [15] 

. d^F^ 3 d^F'' . 

BCD QpBQpCQpD idP^dP^dP^^' 

The following result gives an interpretation of the zero locus of Wilczynski invariants. 

Theorem 2.1 (Grossman [I2j). There is one-to-one correspondence between equivalence 
classes of systems of second order ODEs with vanishing Wilczynski invariants h2.1\) and 
conformal structures of signature (2, 2) with vanishing self-dual Weyl curvature. 

A point p £ M corresponds to an integral curve Lp of the system (jl.ip - a twistor 
line. Moreover, if the Wilczynski invariants vanish then M is endowed with a conformal 
structure and points in T correspond to special null two-surfaces in M called a-surfaces (see 
next section). In Grossman's terminology M acquires a Segre structure from the torsion- 
free system of ODEs |17l EI]. In four dimensions the Segre structure is the splitting of 
the tangent bundle into a product of rank two vector bundles. This is equivalent to the 
existence of a conformal structure on M - see formula (|3.3|) in the next section - so from 
now on we shall use the terminology of conformal geometry without referring to the Segre 
structures. 

In the original Penrose twistor correspondence [23], the four-manifold M is assumed to 
be real analytic, and the correspondence extends to the complexified category. The curves 
Lp in a complex three-fold T corresponding to points in M are globally characterised by 
the type of their normal bundle N{Lp) = 0(1)©0(1), where 0{k) — )■ CF is a holomorphic 
line bundle with Chern class k. In the present context we shall replace this global condition 
on Lp by local-differential conditions (|2.ip . 

3. Twistor Correspondence 

Let [g] = {eg, c : M ^ M+} be a conformal structure on a 4-dimensional manifold M 
consisting of an equivalence class of metrics of signature (2, 2) and let (7 G [^f]. A vector field 
V C TM is called null if giy, V) = 0. The null property of vector fields is invariant under 
rescalings of g, therefore it makes sense to talk about a null vector of a conformal structure. 
Conversely any conformal structure is completely characterised by specifying null vector 
fields. 

Locally there exist real rank two vector bundles S, S' (called spin-bundles) over M equipped 
with parallel symplectic structures e, e' such that 

TM ^ S ® S' (3.3) 

is a canonical bundle isomorphism, and 

g{vi (g) Wi,V2 ®W2) = £{vi,V2)e'{wi,W2) 
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for t'i,f2 G r(S) and 101,11)2 G r(S'). Under this decomposition, any null vector field is of 
the form V = k<^ ir for some k G r(S) and n G r(S'). Also, the Riemann tensor can be 
decomposed as 

Rabcd = i^ABCD^A'B'^C'D' + i^A'B'C'D'^AB^CD 

+<t>ABC'D'£A'B'£CD + (pA'B'CD^AB^C D' 
n 

+t:^{^ac'£bd£A'C'£b'D' — ^ad^bc^A'D'^b'C'), (3.4) 

where ipABCD and ipA'B'C'D' are ASD and SD Weyl spinors which are symmetric in their 
indices and 4>a'B'CD = 4>{A'B'){cd) is the traceless Ricci spinor. 

An a-plane is a two-dimensional plane in TpM spanned by null vectors of the above form 
with TT fixed, and an a-surface is a two-dimensional surface in (" C M such that its tangent 
plane at every point is an a-plane. Examining the Frobenius integrability conditions for 
the existence of a-surfaces leads to the seminal result of Penrose [21]: A maximal, three 
dimensional, family of a-surfaces exists in M iff the Weyl tensor of g is anti-self-dual. 
The anti-self-duality of the Weyl tensor is a conformally invariant property, therefore one 
can talk about ASD conformal structures. Let us assume that the conformal structure 
is real-analytic and therefore can be complexified. Thus (M, [g\) is a holomorphic four- 
manifold with a holomorphic conformal structure. 

Definition 3.1. The twistor space T of a holomorphic conformal structure (M, [g\) with 
anti-self-dual Weyl curvature is the manifold of a-surfaces in M. 

The twistor correspondence puts restrictions on the twistor space T: 

Theorem 3.2 (Penrose i24j). There is a one-to-one correspondence between ASD confor- 
mal structures (M, [g]) and three-dimensional complex manifolds T with the four parameter 
family of sections of ^ with normal bundle 0(1) ® 0{1). 

The additional structures on the twistor space which give Ricci flat metric g C [g] are 

(1) A projection /i : T — > CP , such that the four parameter family of curves above 
are sections of /i. 

(2) Symplectic structure with values in 0(2) on the fibers of ^u. 

The points in the twistor space 7" correspond to totally null self-dual two dimensional 
surfaces (which coincide with the a-surfaces) in M, and points in M correspond to rational 
curves ( in the Ricci flat case these are sections of ^) in T. The relation between 7~ and M 
is best revealed by exploiting the double fibration picture 

M^ F ^T (3.5) 

where the correspondence space 

can be identified with the projective primed spin bundle P(S'). A point in J^ corresponds 
to a point in p G M together with one a-surface containing p. 

The holomorphic curves ^(CP ) where CP = r~^p, p G M, have normal bundle A'' = 
0(1) © 0(1). Two points pi and p2 in M are null separated iff the corresponding curves 
Lpj and Lp^ intersect at one point. 

In the double fibration approach, the twistor space arises as a quotient of J-" by a two- 
dimensional integrable distribution (called the Lax pair) defined by two vector fields Lq, Li 
which at each point of J-" are horizontal lifts to P(S') of vectors spanning an a surface at 
a point in M. The existence of the Lax pair for an arbitrary conformal structure with 
vanishing self-dual Weyl curvature has the following consequence [20^ [T3] . 
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Theorem 3.3. Let Ei,...,E4 be four real vector fields on M and let e^,...,e^ be the 
corresponding dual one-forms. The conformal structure defined by 

g = e^ Qe^ -e^Qe^ 

is ASD if and only if there exists functions /o, /i on M x RP^ depending on X & MF^ such 
that the distribution 

Lo = Ei-XE^ + fo — , Li = E^-\E2 + fi-^ (3.6) 

is Frobenius integrable, that is, [Lo,-Li] = modulo Lq and L\. 

The real (2, 2) ASD conformal structures are obtained by introducing an involution 
r : T — )■ T on the twistor space given by complex conjugation. The fixed points correspond 
to real a-surfaces in M. The involution acts on the twistor lines, thus giving rise to maps 
from CP to CP which swaps the lower and upper hemispheres preserving the real equator. 
The moduli space of r-invariant twistor curves is a real four-manifold. 

4. From ASD Conformal Structures to Systems of ODEs 

In this Section we shall present an algorithm for reconstructing a system (jl.ip of torsion- 
free ODEs from a given ASD conformal structure [g]. Given (M, [g]) we shall construct a 
real projective line parametrised by A G MP of real a-surfaces through a point in p G M. 
This is a curve in M x MP which also depends on the four coordinates of p in M. We shall 
regard this as a parametrised integral curve of a system of two second order ODEs. The 
system (jl.ip will arise by eliminating the coordinates of j* and removing the parametrisation. 
To implement this procedure, consider the Lax pair ()3.6p (see section [3]) of vector fields Lq, 
Li on M X MP . Now proceed as follows 

(1) Find three independent functions {X, y, Z) on F = M x MP which satisfy 

These functions descend to the twistor space T where they provide a local coordinate 
system (X, Y,Z). A point in M corresponds to a curve in T. The pull back of the 
four-parameter family of curves to J- can be parametrised as 

\ ^ [x{\,p),y{\p),z{x,p)), 

where p = {w, z, x, y) are local coordinates on M. 

(2) Use the implicit function theorem to solve the equation X = X[\,p) for A = \{X,p). 
Then solve the relations 

y_V Z-z Y'-^ 7'-^ 

Y-y, z-z, Y - ^^, Z - ^^ 

to express {w,z,x,y) as functions of {Y, Z,Y' , Z'). 

(3) Differentiating once more, and substituting for (w, z, x, y) from above gives a pair 
of second order ODEs (jl.ip . 

As an example of this construction, let us consider the ASD Ricci-fiat metric which, in local 
coordinates may be written as 

g = dwdx + dzdy — Qxxdz — Qyydw +2Qxydwdz (4.1) 

where G^; = dx& etc, and Q = @{w,x,y,z) is a function satisfying the second heavenly 
equation 

Qx^n + Qyz + Qxx&yy " &% = 0. (4.2) 

It follows from the work of Plebaiiski [25] that all (2, 2) ASD Ricci fiat metrics locally arise 
from some solution to this equation. The only non-vanishing part of curvature is given by 
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the ASD Weyl spinor which in the spin-frame adapted to the heavenly equation is a section 
of Sym^(S) — )• M given by 

where the indices A^B, . . . take values 0, 1 and x = (y, —x). The Lax pair for this system 
may be written as 

Lq = dy — X[dw — Qxydy + Qyydx), 

Li = dx + X{dz + Qxxdy-Qxydx)- (4.3) 

A curve Lp C T corresponding to a point p £ M is parametrised by choosing a two- 
dimensional fiber of ^ : T — >• CP and defining {w, z) to be the coordinates of the initial 
point of the curve, and (y, —x) to be the tangent vector to the curve. Thus the pulled 
back curve is A — t- {X = \,Y = y{\,p), Z = Z{X,p)) where the functions {y, Z) admit the 
following expansion |10] 

y = w + Xy-QxX'^ + e,X^ + ... (4.4) 

Z = z - Xx - By x^ - e^ X^ + . . . 

and the higher order terms can be determined by recursion on successive powers of A from 
LA{y) = La{Z) = 0, where A = 0, 1. 

Example 1. Consider a solution to ()4.2p 

e = \y\ 

The corresponding system of ODEs is 

Y" = 0, Z" = -2{Y'f. (4.5) 

The system has a nine-dimensional symmetry group (see formula ()7.ip ). which is also the 
largest symmetry group of a non-trivial ASD conformal structure. The corresponding pp- 
wave metric has constant Weyl curvature with only one non- vanishing component V'oooo = 6. 

Example 2. This example, when analytically continued to the Riemannian signature, 
is relevant in the theory of gravitational instantons. The second heavenly equation ()4.2p 
with 0^ = reduces to a wave equation on a flat (2+l)-dimensional background. Introduce 
t := Qx and perform a Legendre transform 

H{t, y, w) := tx{w, y, t) - @{w, y, x{w, y, t)). 

Then x = Ht, @y = —Hy and ()4.2p yields the wave equation 

Ht^ + Hyy = 0. (4.6) 

together with the metric 

1 ^ T-f 

g = Htt{-dy'^ + dwdt) - -—{dz —dy + Htydwf 

4 ritt I 

= V{jdy'^ + dwdt)-V-^{dz + Af, (4.7) 

where V = Hu and A = Htydw — (Htt/2)dy satisfy the monopole equation *dV = dA where 
* is the Hodge operator on M?'^ with its flat metric. Thus (j4.7p is an analytic continuation 
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of the Gibbons-Hawking metric tl6j . The vanishing of Qz implies that the whole series 
(|4.4|) for y truncates at 2nd order. Prom (|4.4p . we get 

y = w + Xy-X\ (4.8) 

where H = H{w,y,t), from which we can obtain the corresponding path geometry. An 
example with H{w, y, t) = yt^ is 



Y" = X~^Y' - ^/{X~^Y'Y - 2X-1Z', (4.9) 

Z" = {X-^Y' - ^/{X-^Y'Y - 2X-^Z'f/2. 
The potential in the Gibbons-Hawking metric is linear in the flat coordinates on R^'^. 

5. WiLCZYNSKi Invariants 

We shall now demonstrate a new way of deriving the Wilczynski invariants ()2.ip , from the 
double fibration picture ()3.5p . This will provide a converse of the construction presented in 
the last section. Our procedure is analogous to the recursive construction of the Wilczynski 
invariants of a single nth order ODE ^llj . 

Each point in T corresponds to an a-plane in M, and the functions P = [dY/dX, dZ/dX) 
are null coordinates which are mutually orthogonal and thus a^P^ + aiP^ is null for ar- 
bitrary constants a^^ai. Furthermore, by differentiation, the one-form a^dY^ + aidY^^ is 
orthogonal to a^dP^ + aidP^. 

In the derivation below we shall regard {Y ,P ,X) as coordinates on the five dimen- 
sional correspondence space T from the double fibration picture ()3.5p . and define the de- 
generate metric g on this five-dimensional space. We then demand that this quadratic 
form Lie derives up to scale along the total derivative d/dX, and so it gives a conformal 
structure on M. The metric g necessarily takes the form 

g = EABdY^dP^ + ^ABdY^dY^ (5.1) 

where (pAB = 4'(ab) ■ The conformal structure of M is invariant along the fibres of s in the 
fibration ()3.5p . and therefore 

dg ^2 

dX=^' 
for some function Q. Plugging in the expression ()5.ip and comparing coefficients on both 
sides we obtain the equations 

'^(I'AB + Q^eAc = ^^eAB (5.2) 

d(pAB , dF^ 2, ^.oN 

^^+9^^A)C = ^<^AB. (5.3) 

Taking the trace of the ()5.2p (using e to raise and lower indices) 

2Q2 



dF"" .^2 



QpC 

and substituting back into ()5.2p for il we get 

1 dF^ 1 dF^ 

^AB = -2^AC^^ + I^AB-Q^- 

Then, from (|5.3p . we obtain the Wilczynski invariants ()2.ip 



1 d dF^ dF^ IdF^dF^ . 
'2dXdP^^dY^^4dP^dP^^ ^ 
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QpC QpA QpC QpA ^^ 



where here we have used the fact that 

dF^ dl 

dpBdpC - dpCdpB ^ "B 

to rewrite the third term. 

Note here that the expression for the conformal metric (|5.1|) actually gives a metric if 17 = 

i-6, ^pu = which implies 

for some function A. The metric ()5.ip then resembles the heavenly form ()4.ip of the 
Ricci-flat metric. The exact equivalence arises from evaluating ()5.ip at X = 0, where 
Y^ = {w,z),P^ = {-y,x) and A(y^,P^) = -Q{w,z,x,y). 

Example 3. One example is given by 

Y" = Z" = P{Y') (5.5) 

for some arbitrary function /3. The ASD conformal structure on the solution space M is 
type-N and Ricci-flat. Notice here that the system (|4.5p of sub-maximal point symmetry 
is just a special case of this solution. 

The last example was Ricci-flat. One may seek conditions on F for this to be true 
more generally. In other words one seeks to obtain (j4.2p as well as (|4.ip . For this, one may 
proceed as follows: introduce the two- form 

j: = dY AdZ = -EABdY^ A dY^ , 



and then calculate 



and 



%=eABdY^AdP^ 
dX 

= EABidP^ A dP^ + dY^ A dP^), 
"^'^ - EABiSdP^ A dF^ + dY^ A :^dF^) 



dX^ ^^' ' dX 

all modulo dX and where d/dX may be thought of as Lie-derivative. 

In the context of the Nonlinear Graviton Theorem 13.21 the two- form S is the pull back 

of the 0(2)-valued symplectic structure on the fibres of T — > CP to the correspondence 

space. Thus, up to an overall scale, this two-form is given by a quadratic polynomial 

in the base coordinate X. We can therefore examine the consequences of imposing the 

requirement 

d^T, 

„ ~ S modulo dX. 
dX-^ 

The coefficient of dP^ f\dP^ on the left-hand-side must vanish, which imples dP^ /dP^ = 

and we may introduce A as in (j5.4p . Then the vanishing of terms in dY A dP^ reduces to 

a single equation 

.AB 9^A Kad.be g^A d^A ^ 

dY^dpB 2 dP^dP^dP^dpD ''' 



which is precisely (j4.2p in this notation. The remaining condition (j5.3p can now be viewed 
as the initial value problem for A with the initial data given by the solution of the heavenly 
equation. This proves the following 
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Theorem 5.1. Let = Q(w,z,x,y) be a solution to the heavenly equation (|4.2|) which 
gives the Ricci-flat ASD metric (|4.ip . The corresponding system of ODEs with vanishing 
Wilczynski invariants is 

where A|x=Xo = — ©(^(-'^o)) ■^(-'^o)! -^'(-'^o)) — ^'(-'^o)) o,nd the X-dependence of A is de- 



termined by 

The formula (|5.1|) gives a way to reconstruct the conformal structure directly from the 
system of ODEs. An alternative procedure based on the original non-linear graviton con- 
struction of Penrose [23] gives the conformal structure directly from the integral curves of 
the system of ODEs. This procedure will be implemented in the example (j7.2p . 

6. Symmetries of Torsion-Free Path Geometries 

In the twistor approach, we observed there is a one-to-one correspondence between points 
in T and a-surfaces in M. Therefore, any transformation which preserves points in T will 
give rise to a transformation of M which preserves the ASD conformal structure. 

Lemma 6.1. There is a one-to-one correspondence between conformal Killing vectors of 
a (2,2) ASD conformal structure (M, [g]) and point symmetries of the torsion-free system 
of ODEs whose integral curves are twistor lines for (M, [g]). 

Proof. We shall use the double fibration picture ()3.5p . Given a conformal Killing 
vector K of (M, [g]) we can lift it to the vector field K on the correspondence space J-", so 
that [Lo,-ftr] = and [Li,K] = 0, where Lo,Li is the twistor distribution (|3.6p . and the 
commutators vanish modulo a linear combination of Lq and Li. The lift is explicitly given 
hy K = K + Qd\, where Q is a quadratic polynomial in A with coefficients depending on 
coordinates on M. In the spinor notation Q = Vaa'K b'"^ ^ , where vr = (1, —A). 

The space 7" is a quotient if J-' by the distribution spanned by Lq and Li and thus /C = 
s=K-fC is a vector field on 7". Therefore it generates a one-parameter group of transformations 
of T which therefore takes a-surfaces in M to a-surfaces. As K generates diffeomorphisms 
of M and integral curves of the system (jl.ip in T correspond to points in M, then the 
action generated by /C preserves the integral curves of (jl.ip . Thus it is a point symmetry 
on (fLB . 

Conversely, a point symmetry of (jl.ip corresponds to a transformation of M which maps 
a-surfaces to a-surfaces. Therefore it gives a conformal Killing vector on M. 

D 
In the trivial case (Y" = 0, Z" = 0), the path geometry has point symmetry group 5*^(4, M) 
which is isomorphic to the symmetry group of the flat conformal structure 50(3,3). The 
sub-maximal case of a torsion-free system of ODEs with nine-dimensional point symmetry 
algebra (|4.5p which corresponds to a Ricci-flat ASD conformal structure with only one 
constant non-vanishing component of the ASD Weyl tensor. Therefore the 'gap' in (2, 2) 
conformal geometry equals to 6 = 15 — 9. This, by Lemma l6.ll coincides with the gap of 
path geometries in three dimensions, as it is known that the sub-maximal symmetry in the 
projective class V is only eight dimensional. This follows from the general result of Egorov 
|14j that all sub-maximally symmetric projective connections on M" are transitive and have 
symmetry algebra of dimension n^ — 2n + 5. See [19] for an account of the general theory 
of such 'gaps'. 

Example (j4.5p is the unique (up to diffeomorphism) torsion-free path geometry with 
point symmetry algebra of sub-maximal size and we expect systems with 6, 7, 8 symmetries 
to be also comparatively rare. To see explicitly why this occurs, consider the lift of the 
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integral curves of an arbitrary path geometry (|1.1|) to the second jet bundle J^(7', M) which 
is a seven-dimensional manifold with local coordinates given by {X,Y, Z,Y' , Z',Y" , Z"). 
Any point symmetry of p.l|) can locally be described by some vector field x on T which 
we can prolong to a vector field pr^^^x over some open set in J^(7~, K). Then, the functions 

A^ = Y" - F{X, Y, Z, Y', Z'), A^ = Z" - G{X, Y, Z, Y' , Z') 

are constant along pr^^^x- 

Now suppose that the given path geometry admits a Lie point symmetry algebra of di- 
mension five, generated by five vector fields which prolong to an integrable distribution L 
of J^(7', M) on rank at most five. In the generic case, when the rank is five, these vector 
fields will lie in the tangent bundle of some five-dimensional submanifold U C J'^{T,M.). 
Given that the codimension of U is 2, we can construct two functions which A*^, A^ which 
are invariant with respect to the action of L, i.e, given a five-dimensional Lie algebra L 
of vector fields over T, we can find a (non-trivial) path geometry (jl.ip which has L as its 
point symmetry algebra. This statement is true subject to some regularity conditions. For 
example it is sufficient to demand that the first prolongation of the vector fields to the first 
jet bundle 7^(7", M) must not be contained within some four-dimensional submanifold of 
the tangent bundle T{J^{T,R)). 

Lie algebras L of dimension six or greater will generically not give rise to functions 
A'',A^ , and there will be a constraint on finding non-trivial path geometries with point 
symmetry algebras of this size. In particular, we must require that the prolonged algebra 
forms a distribution of rank lower than six. If additionally, we impose the constraint that 
the Wilczynski conditions (|2.ip hold, then this sudden decline of examples will be observed 
sooner (at dimension four rather than six). In the next section, we outline the prolon- 
gation procedure and present some examples of path geometries with 4, 5, 6, 7, 8, 9 point 
symmetries together with some details of the Lie algebraic structure. Where possible, we 
also make use of the twistor approach to say something about the corresponding conformal 
structure on the space of solutions. 

7. Examples 

For a given path geometry in three dimensions (|1.1|) the point symmetries may be found 
by the following well known procedure: 

(1) Let the generator of a point symmetry be 

X = X-idx + XodY + Xidz 

where the Xi ^-re functions of {X, Y, Z) for which we must solve. 

(2) Determine the first and second prolongations by 

^(1) _ ^ _ pA^hc^ (2) _ dv^ _ ^AdX-i 
^^ ~ dX dX ' ^^ ~ dX ^ dX 

(3) The prolongation of the vector field x to the jet bundle J'^{T,M.) is given by 



d sr^ (2) d 



pr 



(^Hx) = x+E^2^7^ + E^i'^ 



B=0 B=0 ^ 



(4) Then we determine point symmetries by finding functions Xi from 

pr(2)(x)(AO)Uo=Ai=o = pr(2)(x)(Ai)Uo=Ai=o = 0. 
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7.1. Special Ricci-Flat Case. Consider the system (|5.5|) with /3 real analytic. The sub- 
maximal torsion-free system (|4.5|) lies in this class so we might expect it to yield more 
examples of systems with a high number of point symmetries. The corresponding ASD 
conformal structure on the moduli space of solutions is given by 

g = dwdx + dzdy — -(3'{y)dw 

and is Ricci-flat, as it corresponds to (|4.ip with 6 = (1/2) / (3{y)dy. 

Under these conditions, it is not difficult to show that the point symmetry algebra contains 
a six-dimensional subalgebra Lq C sl(4, M). It transpires that Lg is solvable and in terms 
of point symmetries we may write it down explicitly 

Lq = spanjei = dx,e2 = 9y, ea = dz, e4 = X^y , 65 = Z9y, eg = Xdx + 2YdY + Zdz}. 

However, as we have seen, there will be some special cases of (j5.5|) for which the point 
symmetry algebra is larger but will contain Lg as a subalgebra. 

Proposition 7.1. Consider a system of two second order ODEs of the form 15. 5|] for some 

function (3 of the form 

00 

fc=0 
If (3 is a quadratic function then the system i5.5\) is diffeomorphic to a trivial one (and the 
symmetry group is 1 5- dimensional) . Otherwise, the symmetry algebra has dimension 

12 - Rank (Mi) - Rank (M2) 
where Mi is a matrix with rows 

(Cfc Ck+i) , k>3 
and M2 is a matrix with rows 

(6 ik-2)^k ik-3)^k-i -(fc + 1)6+1) , k>3. 

Proof. Without loss of generality, let us simplify the problem by making the diffeomor- 
phism 

Z^Z + i^o^' + ^iXY + ^6^' 

so that we obtain the system 

k=3 

which has the same number of point symmetries as the original. 
For a given vector field, 

d d d 

the expressions 

pr(2)(;^)(AO)Uo=A^=o and pr(2)(x)(Ai)Uo=Ai=o 
are real analytic in P^ and P^ with coefficients which are functions of X, Y and Z. For 
X to be a symmetry of the system (4.5), each of these coefficients must vanish separately. 
This leads to the following system of differential equations in x-i, Xo and xi- 

d^Xo d^xi d'^Xo d^Xi d'^X-i d^X-i d^Xo 



dX'^ dX'^ dXdZ dXdY dYdZ dZ^ dZ^ 
^ d'^Xo d'^x-i d'^xo ^d'^X-i d'^Xo d'^x-i 



0, 



dXdY ax2 ' 5y2 dXdY ' dYdZ dXdZ' 
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dXdZ 5X2 ' QYdZ dXdY ' dZ^ dXdZ' 
and, for A; > 3, 

(^ - 2)e.%i + all + (^ - 3)6-i%i - {k + 1)6,,^ - kc,'^ = 0. 

Enforcing the system to be non-trivial (i.e, not all S,k = 0) we obtain 

dx-i ^ 9x0 ^ 9^X-i ^^ 

Then, the most general solution for x is 

X-i = aiX'^ + a2XY + (03 + ai)X + asF + ae 
Xo = 01^^ + a-jX + 02!"^ + 2a^Y + ag 

Xi = aiXZ + a<^X + a2YZ + {az + aiQ)Z + —aiY"^ + (—a-j + aii\Y + ai2 

where, for all fc > 3, 

ika2 +ik+iai = 
and 

&aio + {k- 2)^fca4 + (A; - 3)^^-105 - (A; + 1)^^+107 - {k + l)Cfca3 = 
and the result follows. 

Thus, in the non-trivial case, the point symmetry algebra of (j5.5p has at most dimension 
9. A little more work shows us that systems of dimension 8 are not attainable for this 
example (if the rank of M2 is 1, then the rank of M2 is 2 and if the rank of Mi is 2 then 
the rank of M2 is at least 3). 

n 

Dimension 9. A system with point symmetry algebra of dimension 9 is given by (j4.5p 

Y" = 0, Z" = -2{Y'f. 
A quick check of the expressions for Mi and M2 will show that the system 

also admits a nine-dimensional point symmetry algebra but these systems can be shown to 
be diffeomorphic. The associated Lie algebra in the first case can be written as 

Lg = L6espan{e7,e8,e9} (7.1) 

with 

e7 = Wdv + Zdz , eg = ^Z^dy + Xdz , eg = \x^dx + (]^XY + -^Z^\ dy + \xZdz. 

This has Levi-decomposition 

Lg = Ze K s[(2,R) =span{e2, 63,64,65,67, eg} k span{6i,66 - -67,69}. 
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Dimension 7. A path geometry with point symmetry algebra of dimension 7 is given by 

Y" = 0, Z" = {¥')'' 

for /c > 4. This will be a solvable Lie algebra for any value of k and is obtained by adding 
the vector field 67 = kYdy + Zdz to Lg. 

7.2. ASD Einstein - dimension 8. Consider the system 

Y" = 0, Z" = ^^^ ) ' . (7.2) 

This system has been suggested to us by Boris Doubrov, who has also pointed out its ap- 
pearance in the theory of chains in the homogeneous contact geometry [4j . The Wilczynski 
invariants vanish for ()7.2p . 

The corresponding conformal structure is found by demanding that two neighbouring 
integral curves in T intersect at one point - this is the condition which selects null vectors 
in M. This, in the complexified setting, reduces the normal bundle of the integral curve Lp 
to 0(1) © 0(1). To perform this calculation explicitly, solve ()7.2p for {Y, Z) and demand 
that equations 6Y = and 5Z = for X have a common solution. We apply this procedure 
to the integral curves 

y = w + yX, Z = -+ ^ 

y y\z-xX) 

where 5Y = Yla^^y Sx°-,5Z = Yl,a^a.^ ^^"^i ^"^^ ^'^ ~ {'^■•'^■•^■•V)- Then the discriminant 
condition for two curves to intersect gives the metric 

g = dwdx + dzdy + x dw + [ z -\ \ dy +2\ zx ^ — I dwdy 

\ yj \ yJ 

which is ASD and Einstein, with scalar curvature equal to —24. Both the system of ODEs 
and the metric have eight-dimensional symmetry group S'L(3, M) which acts isometrically 
on M. The homogeneous model for the space of solutions is 5'L(3,M)/GL(2,R), which is a 
(2, 2) real form of the Fubini-Study metric on CP^ = SU{?,)/U{2). 

7.3. Symmetry algebra of dimension 5. The example discussed in the Gibbons-Hawking 
context ()4.9|) has five-dimensional point symmetry algebra which we can write as 

L5 = span{ay , az, A^x + Ydy, -\xdx + Zdz,X^dY + 2Ydz}. 

The algebra L5 is solvable and contains both Bianchi II and Bianchi V as three-dimensional 
subalgebras. 

7.3.1. Sparling-Tod solution. Another way to construct a torsion-free path geometry with a 
given point symmetry algebra is to determine an ASD metric with that conformal symmetry 
algebra and determine the system of ODEs governing twistor lines via the method given in 
Section O The metric of [28] given by 

2 2 

g = dwdx + dzdy — -^ {wdz — zdw) 

[xw + yzy 

is an ASD Ricci-flat metric (14.11) with = — ; — and has five-dimensional conformal 
symmetry algebra. The conformal Killing vectors are 

ki = zduj-xdy,'k2 = -2wdw+xdx-ydy,'k3 = wdw+zdz,'k4, = -yd^+wdz^'k^ = zd^-wdy. 

We should note here that the resulting system of ODEs does not coincide with the previous 
example as the conformal symmetry algebra of the Sparling- Tod solution is not solvable. 
It has an s[(2,M) subalgebra generated by (ki,k2 + k3,k4). Its Levi decomposition can be 
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expressed as the semi-direct product of an s[(2, M) with the 2-dimensional non-Abehan Lie 
algebra. 

The system of ODEs can be read-off directly using Theorem 15.11 and setting A = 
-e(y, Z, Z', -Y'). This yields A = iJ'Z - Y Z')-^ and ([531) gives 

Y" = ^-^ z" = ^^ (7 Q\ 

{Y'Z-YZ'f [Y'Z-YZ'Y' ^ ' 

The integral curves are 

y = Ae^^ + Be-^^ , Z = Ce^^ + De'^^ , 

where {A, B, C, D) are constants of integrations, and 7^ = {AD — BC)~^\/2 . The original 
metric can be recovered from the twistor lines by setting 

w = A + B, z = {C + D), y = -f{B-A), x = -f{C-D). 

7.4. Symmetry algebra of dimension 4. Consider the system 

2 



Y" = o, z" = - z' + V(y')^-i • (7.4) 



This example was found by constructing the most general system of ODEs with Lie point 
symmetry algebra 

L4 = span{dx , dy , dz , Ydx + Xdy} 

and imposing the torsion-free conditions. The algebra L4 above is a particular realisation 
of the abstract algebraU ^4,1 as given by Patera and Winternitz [23]. The integral curves 
of the system ()7.4p are given by 



y = w + Xy, Z = \og{X-x)- ^/y^ - IX + z. 



Following the procedure applied in case of the system (j7.2p we find the ASD conformal 
structure 

g = dxdy + {dw + xdy){dz + yyy"^ — 1 dw). 

This admits a null Killing vector d/dz and thus fits into the classification of ^12j. 

7.5. Symmetry algebra of dimension 3 or less. If the symmetry algebra has dimension 
3, then it necessarily belongs to the Bianchi classification of three-dimensional Lie algebras. 
There are many examples in this case which are analytic continuations of Riemannian 
metrics. See [29] for a discussion of these examples. Examples of an ASD conformal 
class with one or two symmetries can be found in the Gibbons-Hawking class ()4.7p . An 
axisymmetric solution H to the wave equation ()4.6p gives a metric with two Killing vectors. 
Another example with two-dimensional symmetry is the system ()8.5p discussed in the next 
Section. A general solution of ()4.6p with no symmetries gives a metric which only admits one 
Killing vector d/dz. Gravitational instantons of class Dj. are examples of ASD conformal 
structures with no symmetries. 



The Referee has pointed out that there are two more nonequivalent representations of this algebra: 

Lia = span{9y, ~XdY,dz, dx - XZdy} 
and 

L46 = span{29x, 9z, -Y'^dx - ydz, 2Zdx + dy}- 

We have checked that there is no torsion-free system of ODEs with symmetry 1/4;,, and there exists a torsion- 
free system of ODEs with a seven-dimensional symmetry algebra which contains Lia as a sub-algebra. 
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8. FiNSLER STRUCTURES WITH SCALAR FLAG CURVATURE 

For a given n-diinensional manifold T with coordinates X^ {k = l,...,n), a Finsler 
metric is a positive continuous function T -.T ^ [0, oo) such that 

• J^ is smooth on Tr\0 = {(X^ P^') g rr|P / 0}, 

• F{X^,cP^) = c-F(X^P'=) for c> 0, 

• The tensor fij = \-^^ is positive definite for aU {X^,P^) G TT\Q. 

We shall consider the case n = 3 and set X^ = (X, Y, Z) and P^ = X^ = ^^ for some 
parameter t. Finsler geometry generalizes the notion of Riemannian geometry in that the 
norm on each tangent space J-{X , •) is not necessarily induced by a metric tensor. This 
makes these metrics useful in the study of problems involving paths of least time. The 
metric tensor fij allows us to define Finslerian geodesies, which are integral curves of the 
system 

X' + Y^j,X^X'' = 0, i,j,k = l,...,n 
where 

7]fc = 2^^ {flk,X3 + fjl,X'<: - fjk,X') ■ 

It was argued in [5j that, for n > 2, given any system of ODEs with vanishing Wilczynski 
invariants its integral curves arise as the set of unparametrised geodesies of a Finsler function 
of scalar fiag curvature. In this context, the torsion-free path geometries are viewed as 
projective equivalence classes of isotropic sprays on TT- Given a spray 

"'"■'w^-'^w • ' = '•■■■•" 

we can define its Riemann curvature by 

nl fT I V^ I fT I V^ \ J- V^ T"^ P' F™' 

-^kij ~ -"« I '^ jkj -"j I '^ ikj ' ^ im^ jk ^ jm^ ik 

where F*-^ = apjapk ; etc and the Hi form a horizontal distribution determined by S: 



dX^ ' dPi 
Then the spray S is said to be isotropic if its Jacobi endomorphism is given by 

for some function p and covector Tj. It was shown in [6j, that for n > 2, a spray is isotropic 
if and only if it is projectively equivalent to one whose Riemann curvature vanishes (so- 
called R-flat) and in [7], that any R-fiat spray arises as the geodesic spray of some Finsler 
function. More importantly, in [27] it is shown that the geodesic sprays of a Finsler function 
are isotropic if and only if the Finsler function has scalar flag curvature. Here, the flag 
curvature of a Finsler metric is defined in terms of its Riemann curvature by 

i^(x^p^y^^- v\p^r,^mP')v' 



fiP,p)fiv,v)-if{p,v)r 

where y is a section of the tangent bundle transverse to P and the indices are raised/lowered 
with the tensor fij. The flag curvature is scalar if 

K{X\P\V'') = K{X\P''). 

Thus, we could construct all systems of ODEs with vanishing Wilczynski invariants if we 
know how to characterize the Finsler functions of scalar fiag curvature. Although a lot of 
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work has been done in this area, a complete characterization of such functions has not yet 
been achieved. Most success has come with a special type of Finsler function 



F=^a,j[Xk)pipi+h,{X'')P' , ij,k = l,...,n. 

This is known as a Randers metric. If we restrict our attention to Randers metrics of 
constant flag curvature, then we can use the following 

Theorem 8.1 (Bao-Robels-Shen [2].). A Randers metric J- has constant flag curvature if 
and only if the corresponding Zermelo data {h, W) satisfy the following: 

• h is a Riemannian metric with constant sectional curvature. 

• W is a Killing vector or homothety of h. 

Here, the Randers data can be expressed in terms of the Zermelo data as follows [26j : 

_ xhjj + WjWj __Y[i 

aij — 2 'A 

where A = 1 — hijW^W^ and Wi = hijWK This gives rise to a procedure for construct- 
ing torsion-free systems of ODEs from a three-dimensional Riemannian metric of constant 
sectional curvature and a homothety of this metric. The geodesic spray coefficients of such 
systems was worked out in [26|. 

Example 4. Consider the Zermelo data 

d 
h = dX'^ + dY^ + Y^dZ^, W = —-. 

oZ 

The geodesies of the Randers metric associated to this Zermelo data are the integral curves 

of the systems of ODEs with vanishing Wilczynski invariants 



2YZ'^l + {Y'Y + Y^ {{Z'Y - {Y'Y^l) 

y - y(y2 _ 1) (8-5) 

y (1 + {Y'f + {Z'fY'^{{Z'f - {Y'f - 1)) 

+ (y2 _ 1)2 ' 



2Y' z' + yji + {Y'Y + y2 ((^/)2 _ (Y'Y - 1; 
y(y2 _ 1) 

This system has two symmetries -A^ and -^. 

We can also view this correspondence in the other direction i.e, given a system of ODEs 
with vanishing Wilczynski invariants, we can construct a Finsler metric of scalar flag cur- 
vature. 

Example 5. To illustrate this point, let us consider the submaximal system (j4.5p cor- 
responding to ASD Ricci-flat pp waves with constant Weyl curvature. Then, using the 
procedure in [22], this system describes the unparametrised geodesies of the Finsler func- 
tion 

■ (y z xy3 6yy2 \ 
jr = xg\^,2^- 2^— + —. — 
\X X X3 x^ J 

in some open domain of T where fij is positive definite, where Q is any function of two 
variables. In particular, let Q{x,y) = ^Jxy and consider the Lagrangian L = -^J-"^ in 
unparametrised form where X = 1 

L = 2Y'z' - 2X{Y')^ + 6y(y')^. 
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The Euler-Lagrange equations give (|4.5|) . For this example, the flag curvature of J^ vanishes. 

Similarly other systems with vanishing Wilczynski invariants arise from a variational prin- 
ciple induced by a Finsler structure. Thus all these systems fit into the formalism of fl]. 

9. Conclusions and outlook 

We have constructed several examples of systems of two second order ODEs with van- 
ishing Wilczynski invariants. The integral curves of such systems are twistor curves which 
correspond to points in a four-manifold M with ASD conformal structure. The twistor 
curves can also be viewed as unparametrised geodesies of Finsler structures with scalar flag 
curvature. 

The maximally symmetric, non-trivial example is a Ricci-flat ASD pp-wave with 9- 
dimensional group of conformal symmetries, and we have given examples of systems with 
symmetry groups of dimensions between 9 and 5. Some of the examples have a special 
form Y" = 0, Z" = G{X,Y,Z,Y',Z'). This is what Grossman calls 'a weaker form of 
integrability for the second ruling' |,17J. It corresponds to the existence of a two parameter 
family of /3-surfaces (these are null ASD surfaces). This family gives a surface in the twistor 
space T, and it is known [12j any ASD (2,2) conformal structure with a conformal null 
Killing vector admits such structure. 

In the case of trivial ODE with integral curves being straight lines in 7~ = M'^ is a starting 
point for John's integral transform [18j. Given a function i;^ : T — )• K with appropriate 
decay conditions at infinity and a straight line L G T, define a function (p on the space M 
of straight lines 



0(L) 

JL 

The result of John is that the range of this transform is characterised by an ultra-hyperbolic 
wave equation on M, where the wave operator is induced by a flat (2, 2) metric on M = M^. 
If the line L is parametrised by X — )■ {X, Y = w + Xy, Z = z — Xx), then 

4>{w, z, x,y) = / (j){X, w + Xy, z - Xx)dX and 4>wx + (t>zy = 0. 

Jr 

It may be interesting to develop a non-linear version of John's transform applicable to 
non-trivial systems of ODEs with vanishing Wilczynski invariants. Integral curves of (|4.5|) 
could be a good starting point for this construction. 
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